EPJ manuscript No. 

(will be inserted by the editor) 



o 
o 



o : Collective oscillations in quantum rings: a broken symmetry case 



(N 



c3 



c 

o 
o 



> 

O 

o 
o 



i 

c 

o 
o 



X 
S3 



M. Valin-Rodrfguez, A. Puente, and LI. Serra 

Departament de Fi'sica, Universitat de les Illes Balears, E-07071 Palma de Mallorca, Spain 
24 march 2000 

Abstract. We present calculations within density functional theory of the ground state and collective 
electronic oscillations in small two-dimensional quantum rings. No spatial symmetries are imposed to the 
solutions and, as in a recent contribution, a transition to a broken symmetry solution in the intrinsic 
reference frame for an increasingly narrow ring is found. The oscillations are addressed by using real-time 
simulation. Conspicuous effects of the broken symmetry solution on the spectra are pointed out. 

PACS. 73.20.Dx Electron states in low dimensional structures - 78.20.Bh Theory, models and numerical 
simulation 



1 Introduction 

Quantum rings have attracted much interest in recent 
years. They constitute a clear example of the great techno- 
logical advances in the fabrication of electronic nanostruc- 
tures, which also include a rich variety of quantum dots. 
Theoretical studies of rings have addressed the appearance 
of the so called 'persistent current' as a function of an ap- 
plied magnetic field, as well as the electronic structure and 
optical excitations within a fully microscopic approach 
for small rings On the other hand, semiclassical 

methods have been applied to larger rings B, for which 
there exist experimental results in the mesoscopic domain 
H. These semiclassical methods have indeed provided a 
good description of the B-dispersion in the measured far- 
infrared collective excitations for large rings. At an in- 
termediate theoretical level, microcopic density-functional 
calculations have been applied to describe rings containing 
around N sa 10 electrons [HfjJ. An exciting prediction of 
these latter calculations is the transition for increasingly 
narrow rings (approaching the quasi-ld limit) to a state 
with broken rotational symmetry in the intrinsic reference 
frame containing a static spin-density wave (SDW). The 
fabrication of such rings, with a small number of electrons, 
is a very demanding challenge that only recently has been 
addressed by Lorke et a/., who reported measurements of 
the far infrared response of two electron rings J|,[|. 

All calculations of ring excitations for intermediate 
sizes (N > 10), using either semiclassical methods or 
density-functional theory (DFT), have imposed circular 
symmetry to the system and, therefore, have not been 
able to signal the influence of the broken symmetry ground 
state on the spectra. In this paper we relax the symmetry 
condition on the ground and excited states of rings with 
the purpose to explore these effects within DFT. To this 
end we will use real time simulations of the collective den- 
sity oscillations developed by two of us in the context of 



quantum dots |jT^, ^T[. The paper is organized as follows: 
in Sec. II we consider the ring ground-state within our 
approach, focussing on the transition to the symmetry- 
broken ground state; in Sec. Ill we address the real time 
simulation of the charge and spin-density oscillations. Fi- 
nally, the conclusions are drawn in Sec. IV. 



2 DFT description of rings 

In order to model two-dimensional (2d) quantum rings 
of varying widths we have used the model proposed by 
Reimann et al. ^ and consider a ring confining potential 
on the xy plane depending only on the radial distance r, 
such as 



y(ext) ( r 



2 S 



Ro) 



(1) 



In this expression R corresponds to the ring radius and 
we take uj^ — with N p and r s parameters analo- 

gous to those giving for a circular quantum dot the number 
of positive charges and the Wigner-Seitz radius, respec- 
tively. The 5 parameter is introduced in order to model 
the radial thickness of the ring. As delta decreases, the 
ring becomes thinner and thinner, approaching the quasi- 
ld limit. 

We describe the electronic structure within the local- 
spin-density approximation of DFT. The reader is referred 
to Refs. IT2L n3|| for details of the approach in the present 
context. The set of single-particle (sp) orbitals is obtained 
by selfconsistently solving the Kohn-Sham (KS) equations 



_I V 2 + V {ext) {r) + V (H) (r) + V^ xc \v) 



tpi v (r) = £ iv (p iv {r) , 
(2) 



where and vj) xc ^ are respectively the Hartree and 

exchange-correlation potentials. The latter one is obtained 
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from the local energy density as 



y-(xc) 



d£ xc (p, m) 
dp v 



(3) 



In these expressions i] =1,1 labels the spin components 
while the spin densities are p v (t) = J2i IVi^MI 2 - Total 
density and spin magnetization are given by p = p\ + pi 
and m = pi — pi, respectively. Our exchange-correlation 
functional is based on the Tanatar-Ceperley results for the 
2d electron gas [14j with the von Barth-Hedin interpola- 
tion for intermediate polarizations [ p"5[ . 

We emphasize here that our solution to Eq. (||) is not 
assuming a priori any symmetry for the orbitals <p(r). 
The 2d plane is discretized into a uniform grid of points 
in Cartesian coordinates and the Laplacian operator is 
approximated by the corresponding finite differences (we 
have used typically 7 point formulas). The KS solutions 
are then obtained by iteratively applying the imaginary- 
time-step method. This technique is very robust although 
convergence may be rather slow in some cases. The lowest 
energy solution in each case is obtained by performing cal- 
culations with different starting points in order to assure 
that the result does not correspond to a local minimum 
representing an excited state of the system. 



2.1 Circular rings 

Figure 1 shows the density and magnetization for a ring 
with jl6| Rq = 8 , containing N = 10 electrons and for two 
values of the width S. The wide one (S = 2) has perfect 
circular symmetry and its magnetization vanishes. On the 
contrary, the narrow one (5 = 0.08) presents a modulation 
in charge density, a charge density wave (CDW), and spin 
polarization of antiferromagnetic type, with alternate ori- 
entations of spin up and down (a SDW) . This drammatic 
change of behaviour is associated with the increase of the 
radial confinement and agrees with the findings of Ref. 
||. We compare in Fig. 2 the energy of the unrestricted 
solution to Eq. (Q) with that obtained with the contraint 
of circular symmetry, i.e., by solving only the radial KS 
equation as in Ref. Uj. Also shown in Fig. 2 is the ratio 
of maximum magnetization to maximum density and the 
amplitude of the CDW. Panels (b) and (c) show that the 
transition is rather abrupt with the value of 5, although 
the energy (panel (a)) is smoothly changing. Nevertheless, 
the energy trend is very clear and manifests the gain in 
binding energy with the formation of CDWs and SDW's. 

Although some controversy has existed in the quan- 
tum dot community [jl^,[l7|,|l^], the formation of broken- 
symmetry solutions as those displayed by mean-field the- 
ories is a well known phenomenon in nuclear physics Jl9| 
and in atomic physics |^0| . It has been recently discused in 
the context of 2d quantum dots by Yannouleas and Land- 
man 1 18 1 and by Koskinen et al |l^,|l). It has been shown 
that the mean field solution corresponds to the intrinsic 
structure of the system although the full exact solution in 
the laboratory frame preserves the symmetry because of 
an underlying degeneracy. This is convincingly shown by 



comparing internal structure properties of the exact so- 
lution, such as the conditional pair probability, with the 
mean field result [^8|. Another strong evidence of the in- 
ternal symmetry breaking is given by the comparison of 
the low lying states with the roto- vibrational states of the 
structure formed by the localized electrons |l8|,^l). Quan- 
tum rings constitute indeed a very good scenario for this 
behaviour since, as shown in Figs. 1 and 2, they exhibit 
strong breaking of the circular symmetry. A critical com- 
parison of the mean field solution with the exact one for 
two-electron quantum rings will be presented in a future 
publication. 

The broken-symmetry solution is accompanied by a 
strong change in the sp level scheme. While the circular 
structure is characterized by shells, with closings at elec- 
tron numbers of 2, 6, 10, 14 and a small sensitivity on the 
width S, the formation of SDW's and CDWs leads to a 
bunching of the occupied levels and the formation of a rel- 
atively large energy gap at the Fermi level. This behaviour 
is shown in panel (d) of Fig. 2. 



2.2 Elliptic rings 

Deformed rings can be simulated by using an elliptic con- 
tour for Rq(0) of the type 



Ro(0) 



a cos 



1 



■ sin 



(4) 



where a is the major axis of the ellipse, (3 gives the ra- 
tio of minor to major axis and is the polar angle. In 
order to compare with the circular rings of the preced- 
ing subsection we have taken elliptic rings with the same 
perimeter and have considered two deformations, namely 
(3 = 0.75 and (3 = 0.5, corresponding to medium and large 
ring deformation, respectively. 

One of the most striking findings in elliptic rings is 
the appearance of the broken-symmetry solution along the 
ring contour for increasing deformation, at fixed width. 
This happens even for wide rings although the wave struc- 
ture is more clearly marked in the narrow ones (Fig. 3). 
Another behaviour which is apparent in Fig. 3 is the for- 
mation of charge concentrations at the ends of the long 
axis. These are attributed to end states similar to those 
found in finite wires In broad rings, the formation of 
end states appears abruptly after a minimum deformation 
is reached. For instance, at 6 = 2 and (3 = 0.75 no evidence 
is yet found of their formation. 



3 Time dependent spin and charge density 
oscillations 

The analysis of excited states is often made by invoking 
perturbation theory for an external probing field. The fea- 
sibility of this approach usually relies on the symmetries 
of the system, which greatly simplify the treatment. For 
instance, in circular dots one can take density oscillations 
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of the type <5p ?? (r)e zm6 ', i.e., a radial function times a mul- 
tipolar field, and then reduce the equations to the radial 
parts. In the present context the possibility of having a 
broken symmetry ground state forbids this approach. 

A suitable alternative to perturbation theory is to study 
the density oscillations by using real-time simulations. This 
is particularly well adapted to DFT and leads to its time- 
dependent generalization (TDDFT). The simplest version 
of this scheme is the adiabatic-local-density approxima- 
tion, which uses the energy functional of the ground state 
to explore small amplitude oscillations around it. This is 
one of the simplest versions of TDDFT and yet it is quite 
robust and satisfies known exact properties as the energy 
weighted sum rule and the generalized harmonic-potential 
theorem [^2[[2j|. We will use it to calculate the ring spin- 
density oscillations. 

The time-dependent Schrodinger equations can be writ- 
ten 

d 

ig^fir,(r,t) = h v [p,m]<£ iJ7 (r,£) , (5) 

where the KS single-particle Hamiltonian [p, to] is given 
by the square bracket in Eq. (^) , whose stationary solution 
was found in the preceding Section. Here we will perform 
an initial (at t = 0) pertubation of the wave functions 

Pi»j(r,0) =T>(p irl (r) , (6) 

and track the time evolution through Eq. (|J). Technical 
details of the integration method, as well as of the anal- 
ysis of the dipole signals (r)^ which allow to extract the 
excitation spectrum can be found in Rcfs. |Io|.pT|. 

We will distinguish three different types of excitation 
modes, two of them depending on the initial perturbation 
V and another one corresponding to a theoretical model 
in which the electrons are not interacting. Namely: den- 
sity modes associated with an initial rigid displacement of 
the total density; spin modes for which at t = the spin 
densities are rigidly shifted in opposite directions; and free 
modes for which the sp Hamiltonian is kept fixed to the 
stationary densitites, i.e., h v [p,m] « h v [po,mo\. The free 
oscillations, also known as sp modes, are equally excited 
by shifting either the total density or the spin densities. 

3.1 circular rings 

Figure 4 shows the free, density and spin spectra for the 
broad ring with ./V = 10 and 5 — 2. As known from other 
calculations [^Q] the spectrum is roughly divided in two 
regions in the sp model which are then shifted by the inter- 
action in the density and spin channels. The lower region 
contains a single peak while the higher one exhibits an 
important fragmentation. The density response is char- 
acterized by the blue shift from the sp peaks while the 
spin response exhibits a shift to lower energy. The resid- 
ual interaction is also quite effective in enhancing one peak 
above the others in the high energy region. The peak that 
collects the strength is then associated to the more collec- 
tive excitation. These features can be considered standard 
in finite Fermi systems with a repulsive interaction in the 



density channel and an attractive one for the spin channel. 
Indeed, a similar behaviour has been found in quantum 
dots. Figure 4 also shows the corresponding spectra ob- 
tained by using the circularly constrained solution in the 
perturbative random-phase approximation (RPA) |Q. It 
is obvious that both results are equivalent, as one would 
expect since the ground state has circular symmetry. We 
attribute the minor diferences to the quite different tech- 
niques that have been used in the two cases. It is worth to 
notice that the real time method is much more demanding 
computationally and that very fine details of the spectrum 
are harder to describe, such as very low intensity peaks, 
or the resolution of two closely lying excitations. 

The results for the narrow ring (6 — 0.08) are shown in 
Fig 5. Qualitative differences appear with respect to the 
preceding ring and, most important, with respect to the 
circularly constrained solution. They must be attributed 
to the symmetry breaking ground state of this ring. The 
separation of high and low energy region is now clearer, 
with an important gap between both. The most striking 
result comes from the comparison of the low energy peaks 
with the circularly contrained solution. First we notice 
that the full solution in this energy region yields a den- 
sity response which is very close to, but just below, the 
sp peaks and overlaps with some excitations of the spin 
response. In the circularly contrained calculation the den- 
sity response is above the sp one and, furthermore, the 
spin peak lies at an essentially vanishing energy. The dif- 
ferent behaviour is indicating that the symmetry breaking 
is giving stability to the ring, which would otherwise be- 
come unstable against spin oscillations. This is in agree- 
ment with the appearance of an energy gap in the sp en- 
ergies associated to the symmetry breaking, mentioned in 
Sec. II. 

The high energy region is associated to a locally radial 
oscillation in the ring, and thus similar to the oscillation 
in a one dimensional oscillator. This is confirmed by the 
similarity of the energy with the radial curvature at the 
minimum given by that, by the well known gen- 

eralized Kohn's theorem, is the only allowed energy for 
rigid density oscillations of dipole type in parabolic con- 
finement p4| . On the other hand, the low energy region is 
associated with oscillations locally tangential to the ring. 
Within this interpretation, the fact that the density and 
spin tangential modes are close in energy and essentially 
overlap with the sp excitations can be understood as a sig- 
nature of a 1-dimensional Luttinger liquid p5[ , for which 
the lowest excitations are not the sp but the collective 
ones. 

In the symmetry broken ground state there is also a 
possibility to excite pure tangential modes by means of 
a twist of the two spin densities in opposite directions. 
These exotic modes were proposed recently in the context 
of quantum dots (h}. It is obvious that pure tangential 
modes can not exist in circularly symmetric rings since a 
twist of spin densities does not change the energy. How- 
ever, in a symmetry broken state a restoring force appears 
because of the distortion of the spin-density wave. The as- 
sociated oscillation energy is easily obtained by analyzing 
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the time dependence of the orbital currents {Y^iH a i ) 
that appear after an initial spin twist. 

Figure 6 shows the spin twist spectrum for the narrow 
ring with 10 electrons. We observe that the peaks are in 
the same region as the low energy dipole excitations, which 
is proving that the latter ones are indeed of tangential 
nature. Both free and interacting spectra are characterized 
by a rather regular energy spacing, after an initial gap. We 
attribute this to the parabolicity of the potential in the 
tangential direction for points close to the local minima. 

3.2 Elliptic rings 

The systematics with deformation is very similar to that 
discussed for the circular case, with the additional result 
that both dipole energy regions are now more fragmented. 
Figure 7 shows the corresponding spectra for a strongly 
deformed ring (3 — 0.5 with thicknesses given by S = 2 
and S = 0.08. This fragmentation is more regular in the 
wide rings, with a doubling of high and low energy peaks 
by the same energy separation. In the narrow ring the 
relative fragmentation is lower with respect to the separa- 
tion of high and low energy regions which is in agreement 
with the previously mentioned transverse and longitudinal 
character of both regions. 



Kohn's theorem. The separation between high and low 
energy regions incresases as the ring width decreases. 

Pure tangential modes in rings with broken symmetry, 
associated to a twist of the spin densities have been dis- 
cussed and it has been shown that their energy is close 
to the low dipole peaks. This supports the interpretation 
of this dipole excitations. Besides, we have shown that 
the more exotic spin-twist modes are regularly spaced in 
energy because of a high degree of harmonicity of the po- 
tential. Finally, the dipole oscillations in elliptic rings have 
been obtained and their fragmentation depending of the 
ring width has been shown. 

Further work is necessary in order to analyze in detail 
the convergence to the strict Id limit, for very narrow 
rings, and how the peculiarities of a Luttinger liquid such 
as power-law (as opposite to long range) order in the SDW 
correlation function and spin-charge separation p5[ | may 
manifest in a finite size system. 
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4 Conclusions 

The ground state and oscillations of quantum rings have 
been discussed within density-functional theory relaxing 
the constraint of circular symmetry of the electronic den- 
sity. A transition to a broken symmetry ground state with 
static charge and spin density waves has been found for 
increasingly narrow rings, in agreement with previous re- 
sults. The energy gain with respect to the circularly con- 
traint solution, the amplitude of charge and spin-density 
wave as well as the sp level scheme have been discussed as 
a function of ring width. The transition to a wave structure 
along the ring contour has also been obtained in elliptic 
rings for increasing deformations. In this case, there also 
appear charge concentrations at the long axis ends indicat- 
ing the formation of end states similar to those obtained 
in finite quantum wires. 

Signatures of the broken symmetry ground state on the 
collective oscillation energies have been searched for in the 
dipole spin and density oscillations. We have shown that 
the broken symmetry is accompanied by an overlap of the 
density and spin response peaks in the low energy region 
corresponding to tangential oscillations. At the same time, 
the low lying sp excitations are at a slightly higher energy, 
which is indicating a quasi Id behaviour similar to that of 
Luttinger liquids. The energy of these excitations is close 
to the gap associated to the broken symmetry. It has also 
been shown that the broken symmetry gives stability to 
the ring, which would otherwise become unstable against 
spin dipole oscillations. The oscillations in the radial di- 
rection are located at higher energy, close to the radial 
curvature of the parabola as indicated by the generalized 
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Fig. 1. Results for a ring with N = 10 electrons, Ro — 8 and 
parameters N p = 10, r s = 1.51 (see text). From left to right 
the density for a ring width 8 = 2, the density for 8 = 0.08 and 
the magnetization for 8 = 0.08 are shown, respectively. 
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Fig. 2. Evolution with 8 of several properties for the same ring 
of Fig. 1. Panel (a) shows {E C ~E) / E in percentage, where E c is 
the energy of the circularly contrained solution while E is that 
of the full solution to Eq. (2) . Panel (b) shows the evolution of 
the ratio of maximum magnetization over maximum density. 
Panel (c) displays the amplitude of the charge-density wave 
and panel (d) shows the evolution of the energy level scheme 
referred to the lowest sp orbital. In this latter panel, occupied 
levels are plotted with solid lines while the first unoccupied one 
is shown with a dashed line. 



Fig. 3. Density and magnetization for elliptic rings with 10 
electrons and 8 = 2. For /3 = 0.75 only the density is shown, 
since the magnetization vanishes, while for f3 = 0.5 both den- 
sity and magnetization are displayed. 
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Fig. 4. Dipole oscillation spectra in the free, density and spin 
channels for the circular ring with N = 10 and 8 = 2. The 
left panels correspond to the unrestricted solution while the 
right panels show the circularly constrained one. The arrow 
shows the value u>o/V8, that gives the curvature of the radial free and interacting channels, 
parabola at the minimum. 



Fig. 6. Spin twist spectra for the narrow circular ring in the 
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Fig. 5. Same as Fig. 4 but for 8 = 0.08. 



Fig. 7. Same as Fig. 4 for a deformed ring with f3 = 0.5. Left 
panels correspond to 8 = 2 and right ones to 8 = 0.08. 
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